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Abstract. We present the results of two-dimensional magnetohydrodynamical nu- 
merical simulations of relativistic magnetic reconnection, with particular emphasis on 
the dynamics of Petschek-type configurations with high Lundquist numbers, S ~ 10^ - 
10*. The numerical scheme adopted, allowing for unprecedented accuracy for this type 
of calculations, is based on high order finite volume and discontinuous Galerkin meth- 
ods as recently proposed by Du mbser & Zanotti ( 2009 ). The possibility of producing 
high Lorentz factors is discussed, by studying the effects produced on the dynamics by 
diff'erent magnetization and resistivity regimes. We show that Lorentz factors close to 
~ 4 can be produced for a plasma magnetization parameter cr,„ = 20. Moreover, we 
find that the Sweet-Parker layers are unstable, generating secondary magnetic islands, 
but only for S > S c ~ 10*, much larger than what is reported in the Newtonian regime. 



1. Introduction 



Relativistic magnetic reconnection is a high-energy process converting magnetic energy 
into heat and plasma kinetic energy over short timescales. It is supposed to play a 
fundamental role in the magnetospheres of pulsars ( [Uzdensky 2003 Gruzinov 2005| ); at 
the termination shock of a relativistic striped pulsar wind (Petri & Lyubarsky 2007); in 
soft gamma-ray repeaters (,Lyutikov.2003^ 2006^; in gamma-ray burst jets (,Drenkhahn 



& Spruit||2002HBarkov & Komissarov||2010HMcKinney & Uzdenskyl|20T0 
et al.'201 1 1; in accretion disc coronae ( diMatteo 1998 Schopper et al.|1998" 
eral. |2004 . 



RezzoUa 



Jaroschek 



In spite of the initial optimistic expectations that relativistic magnetic reconnection 
could provide very fast reconnection rates, evidence has emerged over the years that the 
relativistic Petschek reconnection should not be considered as a mechanism for the di- 
rect conversion of the magnetic energy into the plasma energy and the reconnection 
rate would be at most 0. 1 the speed of light, contrary to what originally suggested by 
Blackman & Field ( 1994[ l. Fundamental progresses in the numerical modeling of rela- 
tivistic magnetic reconnection have been recently obtained by Watanabe & Yokoyama 
(2006); Zenitani et al.| ( [2009 b a), However, two major numerical limitations still pre- 
vent realistic astrophysical applications of numerical schemes specifically devoted to 
relativistic magnetic reconnection. The first limitation is due to the difficulty in reach- 
ing sufficiently high magnetization parameters cr,„, while the second limitation is due 
to the difficulty in treating physical systems with very high Lundquist numbers S . 
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In this work, by adopting the innovative numerical method presented in Dumb- 



ser & Zanotti (12009), we show the results of numerical simulations in the very high 
Lundquist numbers regime, S ~ 10^-10**, showing that that Sweet-Parker current 
sheets are unstable to super- Alfvenically fast formation of plasmoid chains, but only 
for 5 >Sr-' 10^ 



2. Physical set up and numerical approach 

The initial model that we have considered is built on Harris model, as reported by 



Kirk & Skjaeraasen (2003 1, and it reproduces a current sheet configuration in the x - y 
plane. Gas pressure and density are given by p - po + cr^pQ[pocosh^{2x)]'^ , p - 
Pq + CTmPolpo cosh^(2x)]"', where po and po are the constant values outside the current 
sheet, whose thickness is 5 = 1. The magnetic field changes orientation across the 
current sheet according to By = Botanh(2x), where the value of Bq is given in terms 
of the magnetization parameter o",„ = B'^^KlpoF'^). All over the grid there is a small 
background uniform resistivity ijh, except for a circle of radius = 0.8, defining a 
region of anomalous resistivity of amplitude 7/,o = 1 0. The resistivity can be written as 

T]t + 7]io [2(r/r,,)^ - 3{r/r,jf + l] for r < r^, 
m for r > rr^, 

where r = yjx^ + y^. The velocity field is initially zero, while the electric field is given 
by fj. = 7]{dByldx). We have considered the case with /^o - 1» Po - 1- The Lundquist 
number for every model is 5 = vaL/t/^, where L is the length of the initial current sheet, 
while v\ = B^/{hp + B^) is the relativistic Alfven velocity. 

A well known and challenging feature of the relativistic resistive magnetohydro- 
dynamics equations is that the source terms in the three equations for the evolution 
of the electric field become stiff in the limit of high conductivity. To cope with this 
difficulty, we have applied the strategy described by |Dumbser & Z anotti (2009), who 
used the so called high order PnPm methods, which combine high order finite volume 
methods and discontinuous Galerkin finite element schemes in a more general frame- 
work (Dumbser et al. 2008b a). The numerical grid consists of an unstructured mesh 
composed of triangles which are clustered along the current sheet. The grid extension 
is given by [-50, 50] x [-150, 150]. We have used periodic boundary conditions at 
and jmax, while zeroth order extrapolation is applied at x^m and Xmax- 



3. Results 

The dissipated magnetic energy, triggered by the anomalous resistivity, produces an 
increase of both the thermal and kinetic energy. The latter results in the acceleration of 
the plasmoid along the direction of the magnetic field, which is more efficient for higher 
magnetizations. However, as the region around the anomalous resistivity becomes more 
and more rarefied, the conversion of magnetic energy into thermal energy becomes 
more efficient than the conversion into kinetic energy and the Lorentz factor reaches a 
saturation. This effect is shown in Fig. [T] 

Fig. [2j on the other hand, shows the development of an tearing-like instability 
for very high Lundquist numbers configurations. No sign of instability is visible in 
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Figure 1 . Time evolution of the Lorentz factor for models with increasing mag- 
netizations (Tm and S ~ 10^. 




Figure 2. Generation of plasmoid chains in very high Lundquist numbers config- 
urations. Left and right panels report, respectively, the color map of the rest mass 
density at time t = 60 for two models having S = 5.36 x lO'' and 5 = 2.68 x 10^ 
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simulations with Lundquist numbers as large as 5 ~ 10^ (left panel), while a chain of 
magnetic islands is produced for 5 > 5^ ~ 10*^ (right panel). Our results, combined 
with those by Samtaney et al. (2009 1, indicate that, in the transition to the relativistic 
regime, the critical Lundquist number increases from S c ^ 10** to > 10*^. Such 
a conclusion may have deep implications for high Lundquist number reconnection in 
realistic astrophysical conditions, and it requires further investigations. 



4. Conclusions 



By adopting high order discontinuous Galerkin methods as proposed by |Dumbser~& 



Zanotti (2009), we have found that Lorentz factors up to ~ 4 can be obtained for plasma 



parameters cTm up to 20 in systems undergoing relativistic magnetic reconnection with 
Lundquist number S ~ 10^. When S is larger than a critical value S c ~ 10**, the 
Sweet-Parker layer becomes unstable, generating a chain of secondary magnetic is- 
lands ([Zanotti & Dumbser|201 
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